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Polar optical phonons are studied in the framework of the dielectric continuum approach for a 
prototypical quantum-dot/quantum- well (QD/QW) heterostructure, including the derivation of the 
electron-phonon interaction Hamiltonian and a discussion of the effects of this interaction on the 
electronic energy levels. The particular example of the CdS/HgS QD/QW is addressed and the 
system is modelled according to the spherical geometry, considering a core sphere of material "1" 
surrounded by a spherically concentric layer of material "2" , while the whole structure is embedded 
in a host matrix assumed as an infinite dielectric medium. The strength of the electron-LO phonon 
coupling is discussed in details and the polaronic corrections to both ground state and excited state 
electron energy levels are calculated. Interesting results concerning the dependence of polaronic 
corrections with the QD/QW structure size are analyzed. 
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I. INTRODUCTION 



Quantum-dot/quantum-well (QD/QW) semiconduc- 
tor heterostructures based on CdSe, ZnS, CdS, HgS, 
etc., have been synthesized by means of colloidal solu- 
tion chemistry in the latter years. They display rather 
interesting size-dependent optical propertiesi2ii^i& and 
appear to have promising perspectives in the streaming 
progress of nanostructured device technology. In the case 
of QD/QW heterostructures a few monolayers of ma- 
terial "2" (say, /3-HgS) are epitaxially grown on a core 
structure of material "1" (say, CdS) and, after that, the 
system is sometimes capped with an additional layer of 
material "1". The whole system is embedded in a host 
material, which may be a silicate matrix or an organic 
compound. We should remark that, for QD/QW het- 
erostructures based on CdS//3-HgS, a true quantum- well 
is created within the QD because the bulk /3-HgS energy 
gap is lower than the CdS gap. Highly confined carriers in 
the /3-HgS layer were actually detected by means of differ- 
ent kinds of optical experiments. Besides the electronic 
and optical properties, the vibrational spectra are rather 
important in the discussion of the physical characteristics 
of these heterostructures. Usually the so-called dielectric 
continuum approach (DCA) has been applied for the the- 
oretical determination of TO and LO bulk-like phonons, 
while, in the case of small-size crystal heterostructures, 
we must also consider the SO phonons associated to 
the system interfaces^. Most of the theoretical studies 
considered the spherical geometr y 7 ! 8 : 9 : 10 ! 11 , but in some 
cases spheroidal QD's were investigatecU2iiii4. Theo- 
rctical ly more reliable approaches, as those reported in 
Refs. |15llfill7| . involve rather complicated calculations. 
However, it must be realized that the DCA, being an es- 
sentially electrodynamical approach and relatively simple 
in its theoretical formalism, is still able to reveal many 
important physical features of the system as has been 



proved by all previous works on the subject. 

In this paper we focus on the confined internal LO 
phonons, and also on the electron-LO phonon interac- 
tions for a QD/QW heterostructure. The corresponding- 
analysis for SO phonons was already accomplished in Rcf. 
[pc| for the same kind of model-system. By applying the 
methods of the DCA we derive the electric potential for 
confined internal LO phonons and, using this result, the 
electron-LO phonon interaction Hamiltonian is deduced. 
The electric potential allows to estimate the strength of 
the interaction, but the interaction effects may also be 
estimated by an analysis of the energy shifts of the unper- 
turbed levels (the so-called polaronic effect). We discuss 
that this kind of effects is able to provide a better esti- 
mation of the electron-phonon interaction strength. Po- 
laronic effects have been studied by different authors for 
systems more or less similar to the one we are discussing 
in the present paper. Within the adiabatic approxima- 
tion, and using second order perturbation theory, the po- 
laronic corrections were studied for nanocrystallites in 
Refs. |ToKl20j . and also in Refs. |TlTl2lI^l2l| by apply- 
ing other methods. It is worth to emphasize that nanos- 
tructured QD's, considering systems beyond the single 
QD sphere, involve more complicated vibrational prop- 
erties and need to be further studied. Henceforth, in 
the present work we aim to show the main effects of the 
electron-LO phonon interactions and their influence on 
the energy levels of the QD/QW heterostructure by ap- 
plying a relatively simple model. We show calculations 
involving both the ground state energy level and the ex- 
cited state energy levels of the QD / QW heterostructure 
and comparisons are made with other previous works on 
the subject. 

The paper is organized as follows. In Sec. II we give a 
brief summary of the fundamental results of the DCA and 
apply them to our model-system, thus deriving the inter- 
nal LO phonon electric potential. In Sec. Ill we analyze 
the wave functions and energy levels in the framework of 
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the effective mass approximation and give a perturbative 
estimation of polaronic energy shifts. Finally, Sec. IV 
presents a detailed discussion of the main results of the 
work. 



II. ELECTRIC POTENTIAL FOR THE 
INTERNAL LO PHONONS 

The DCA is based on the well-known Born-Huang 
equations, which have been extensively described in the 
Iiterature£i£i2ii2iiii and leads to the equation: 

e(w)VV = 0, (1) 

where (p is the quasistatic electric potential and 
s(lo) is given by the standard expression: e(y>) = 
Sao (lu 2 —uif) I (lu 2 — uij,). In the above equation the 
harmonic time dependence <p(r, t) ~ exp(— iujt) is ap- 
plied. We further assume the Lyddane-Sachs- Teller re- 
lation: lu^/lu^ — Sq/Eoo, while £o (£00) is the static 
(high frequency) dielectric constant of the medium and 
u>l, wt are the bulk longitudinal and transverse phonon 
frequencies at the semiconductor T point. The internal 
LO phonons satisfy Eq. (JIJ with e(lu) = and, therefore, 
lu = lul . Obviously, the electric potential does not satisfy 
the Laplace equation. Thus, the potential must be ob- 
tained as an expansion in an appropriate system of base 
functions fulfilling the fundamental symmetry conditions 
of the structure. In our case, the structure is spherically 
symmetric, and the spherical coordinates (r, 0, <fi) are the 
more convenient. Of course, the electric potential must 
also satisfy the usual boundary conditions: p should be 
continuous at the interface between two different media 
and also should satisfy the condition 
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i.e., the continuity of the normal component of the elec- 
tric displacement D. Let us also recall the rather useful 
relation 

Vc/?= — W— w, — = , (3 

Eoo - £{lu) V £* £* £oo £o 

where p is the reduced mass density addressed to the 
involved ion couple and u is the mechanical relative dis- 
placement vector with units of length. For the internal 
LO phonons Eq. J3J| is applied taking e(uS) = 0. 

The geometry of the spherical QD/QW heterostruc- 
ture is assumed in the following form: in the interval 
< r < a we have material "1" (a polar semiconductor 
like CdS), in the interval a < r < b we have material "2" 
(a different polar semiconductor like HgS), and finally 
for r > b we shall assume an infinite dielectric medium 
with a fixed dielectric constant Ed- Using spherical co- 
ordinates the electric potential is expanded as follows: 
for r < a we apply the complete set of orthogonal func- 
tions ji(qr)Yi m (0, </>), where ji(x) with I = 0, 1, 2, • • • are 



spherical Bessel functions and Yi m (8, (f>) are the spher- 
ical harmonics with \m\ < I; for a < r < b we shall 
apply the complete set of orthogonal functions defined 
by [ni(qa)ji(qr) - ji(qa)ni(qr)}Yi m (0, <fi), where rii(x) are 
spherical Bessel functions of the second kind. Let us 
remark that in the present model-system two different 
semiconductor materials are involved and two different 
types of phonons should be taken into account. In the 
central spherical core (r < a) we have semiconductor "1" , 
oscillating with the frequency loil, in the spherical layer 
(a < r < b) we have semiconductor "2" oscillating with 
the frequency t02L- For the same reason we must con- 
sider two different electric potentials associated to these 
phonons: electric potential due to phonons in the 

core sphere, and electric potential p^ 2 ' due to phonons in 
the concentric spherical layer. There is not overlapping 
between these electric potentials: p^ is confined to the 
core sphere and must be zero for r > a, while 

pW i s 

confined to the spherical shell and must be zero for r < a 
and r > b. The latter properties of the electric poten- 
tial are a clear consequence of the boundary conditions 
(|2*|l when applied to the two spherical surfaces at r = a 
and r — b respectively. On the other hand, the equation 
jl (qa) — determines the possible values of q involved in 
the expansion of the electric potential p^\ For the elec- 
tric potential p^ the possible values of q present in the 
corresponding expansion are determined by the equation 
ni(qa)ji(qb) = ji(qa)ni(qb). 

It is useful to represent the electric potential in the 
following form: p^{r, 0, (f>) — ^\ , '\r)Yi m (0, 0), where 

<&P(r) (i — 1, 2) is the radial part of the potential, which 
is also dependent on the various values of q defining the 
possible phonon modes. It is then clear that for r < a 

*?\r) = 4%)3i(ir), (4) 

and 

®i 2) (r) = A[2(g)[n,(ga)ji(9T)-J«Mni(«r)], (5) 

for a < r < b. In Eqs. (QJ and JSJ (q) are normaliza- 
tion constants (to be determined below) and the expan- 
sion involves summations over q, /, and m. In order to 
proceed with the normalization of the electric potential 
we should introduce the phonon annihilation and cre- 
ation operators ai m (q) and ctj (q) obeying bosonic com- 
mutation relations. Correspondingly, the classical dis- 
placement vector u is transformed into the operator u 
by means of 

fiW =4 i} V/W(r)a Im (g) (6) 

where 

/£» = Mqr)Y lm (e,4,), (7) 

/£» = Mqa)ji(qr)-Mqa)ni(qr)]Y lm (9 t 4i).(8) 

We should realize that the q's present in potentials of 
type "1" or "2" are different (determined by different 
transcendental equations) . 
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In order to determine the constants Uq a well-known 
procedure should be applied, which has been described 
elsewhere2*i2iiii&. We transform the classical kinetic 
energy due to the vibrations into a quantum operator, 
which is then taken as the LO phonon Hamiltonian H p h- 
Requiring this Hamiltonian to have its standard har- 



monic form 
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hujiL^a + 1/2) we can deter- 
mine Uq' . The calculations are straightforward and we 
briefly show them in the Appendix. We have found that 
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where e is the absolute value of the electron charge. In 
Eq. (|16|) we have included a summation over all the pos- 
sible eigenfunctions of this model-problem. In some ex- 
pressions we are using simultaneously the notations q or 
x = qa in order to denote the roots of the corresponding 
transcendental equations. Whenever the notation xi is 
applied the roots of the corresponding modes for i = 1, 2 
are implied. This should not lead to confusions. 



III. ELECTRON ENERGY LEVELS AND 
POLARONIC CORRECTIONS 

Applying the effective mass approximation for the elec- 
tronic states, we model the confining potential for the 
electron as 



where 

Mi{x) = 7 3 [ji(x)n l+1 (-yx) - ni(x)ji +1 (-yx)} 2 

- \ji(x)m+i(x) - ni(x)ji +1 (x)} 2 , (11) 

and x = qa are the roots of the corresponding transcen- 
dental equation: ji(x) — for phonons of type "1" and 
ni(x)ji(^x) = ji(x)ni(jx) for phonons of type "2". We 
have also introduced the parameter 7 = b/a. 

By the above described method we finally obtain the 
normalization constants A^lq) (notice that actually 
they do not depend on m due to the azimuthal sym- 
metry of our model) introduced in Eqs. (0} and (JSJ) . In 
this step Eq. © must be applied and we obtain 



U(r) 



-U 



for < r < a 
for a < r < b 
for r > b 



(17) 



where Uq is the band-offset between the conduction bands 
of semiconductors "1" and "2". We are here assuming 
that semiconductor "2" has a smaller energy gap and, 
therefore, in the shell layer a < r < b the conduction 
electrons are actually subjected to a QW. The zero en- 
ergy level was chosen at the bottom of the semiconductor 
"1" conduction band. 

The solution of the Schrodinger equation H^!(r) = 
E^(r) is a straightforward problem and we are led to 
the wave function 
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In Eqs. U21 and 1/e* = l/e ioo - l/e i0 . We have 
thus obtained an explicit analytic expression for the elec- 
tric potential operator ip^ corresponding to both types 
of internal LO phonons involved in our model-system: 



^l(r, 8, 0) 
<P?L(r, 9, 0) 



\Af\xi) [f$l(r)a im (x. 2 ) + ft.cjlS) 



where 'h.c.' stands for 'hermitian Conjugate'. Let us 
once again emphasize that the electric potential operator 
with i = 1 is defined as the zero operator for r > a, while 
the operator with i = 2 is the zero operator for r < a 
and r > b. The electric potential operators defined in 
Eas. l|14|l and (|15|) are hermitian and allow us to define 
the electron-LO phonon interaction Hamiltonians, which 
are given by 
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(16) 



*(r, 9, <j>)=i ! {r)Y lm {6, 



where 



ip(r) = C ■ 



ji(k r) 



ni (kb)ji (kr)-ji (kb)ni (kr) 
™i ( fc& ) 3l (ka)-ji (feb)ni (ka) 



In the above equation we have defined 



k = 
k = 



l2fi 2 (Uo + E) 



(18) 

for r < a 
for a < r < b 
(19) 

(20) 
(21) 



qlm 



where with i = 1, 2 denote the electron effective mass 
for each semiconductor. In this context the symbol po 
stands for the free electron mass. The energy levels, de- 
noted by E, are all discrete and could be positive or nega- 
tive. For negative energy levels (always larger than — Uq) 
we should take ko = in, where k is again given by H20(l 
taking \E\ in the place of E. The normalization constant 
C can be calculated requiring J \^\ 2 d 3 r = 1. 

Notice that the wave-function described in (|19fl is con- 
tinuous at r = a and equal to zero at r = b. The addi- 
tional boundary condition 



4 



1 




1 


~dip2~ 




Ml 


dr 


a A*2 


dr 


a 



leads to the discrete energy levels, given by the roots of 
the transcendental equation 

fco^2 j/Cjgpo) _ ni(kb)j{(ka) - ji(kb)n[(ka) 
kni ji(k a) ni(kb)ji(ka) - ji(kb)ni(ka) 

In Eq. (J23J the prime in the spherical Bessel functions 
denotes the first derivative with respect to the function's 
argument. 

It is interesting to analyze the effects of the electron- 
LO phonon interaction on the electron energy levels. For 
the considered semiconductors this interaction is, in gen- 
eral, weak and should lead to a relatively small shift to- 
wards the lower energies. Taking into account just the 
interaction with confined internal LO phonons we realize 
that, for the studied QD/QW heterostructure, two dif- 
ferent electron-LO phonon interaction Hamiltonians are 
involved, as predicted by Ea . (fltjfl . For LO phonons of 
type "1" ("2") only the region r < a (a < r < b) provides 
a contribution. In general, the SO phonons (analyzed in 
Ref. ^8|) a l so contribute to the energy shift for all the 
electronic states with I > 1. For electron states with 
I = the SO phonons do not contribute to the energy 
shift. Other important issue is the degeneracy of electron 
states with I > 1. This degeneracy (of order 21 + 1) in 
general could be lifted if the interaction with LO and SO 
phonons is taken into account. In the present work spin 
effects arc not considered. 

The general theory for the calculation of the polaronic 
corrections in systems similar to the one studied in this 
work has been described elsewhereii^i^iSSiSiiS^ and need 
not be given here in detailed form. Applying the pertur- 
bative approach (and assuming that the system behaves 
adiabatically) we are led to the expression 

^\{p,l\H?l ph \p, 0) 2 

In Ea. (|24|l we are assuming T — K, so all the LO 
phonons involved in the polaronic corrections are vir- 
tual. Label "p" is used to describe the electronic state 
we are analyzing, while label "p"' denotes all the in- 
termediate electronic states in the second order pertur- 
bation theory formula. AE p gives the energy shift of 
level "p" due to the electron-phonon coupling. The sym- 
bol Ep corresponds to unperturbed electron energies and 
\p, 1) (|p, 0)) is a ket- vector notation for the electron-one 
phonon (electron-zero phonon) states. In our concrete 
calculations the symbol "p" involves three quantum num- 
bers (I, m, v) but the energy levels depend just on '7, v" , 
where 'V denotes the roots of the corresponding tran- 
scendental equation (Eq.J^BJ)) for a fixed value of I. If 



we are renormalizing the energy of an electron state with 
/ = the calculations become easier because the cor- 
responding level is not degenerate and the contribution 
from SO phonons can be ignored. States with I > 1 are 
obviously more complicated and require quite involved 
calculations. Another important issue present in Eq.J23J) 
is the summation over all possible phonon states. How- 
ever, it should be kept in mind that the possible phonon 
modes are, in general, bounded, i.e. the number of pos- 
sible modes cannot be larger than the system volume di- 
vided by the unit cell volume. In the practical case, just 
the first phonon modes provide a significant contribution 
to the energy shift. 

IV. DISCUSSION OF THE RESULTS 

Let us discuss the essential results obtained in the fore- 
going sections. For the numerical computations we con- 
sidered that material "1" is CdS and material "2" is (3- 
HgS. The values of the corresponding physical parame- 
ters are given in Table I. 



Table I: Semiconductor physical parameters. 
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9.1 


233 


300 


0.2 


2.5 


HgS 


11.36 
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250 
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We must emphasize that, within the adopted model, 
the properties of the dielectric host medium for r > b 
are immaterial. Both internal LO phonons and electrons 
are defined for r < b. SO phonons, however, do depend 
on the properties of the host medium, but they are not 
considered here. 

One important issue which deserves a close analysis 
is the strength of the electric potential associated to 
the two different types of internal LO phonons of our 
system. A first estimation can be made by a study 
of the radial part of the electric potentials as given in 
Eqs. Q and JHJ|. In Fig. 1 we are plotting <&W(r) 
(in units of $ = (4?iwiL/7rae*) 1 / 2 ) for I = 0, setting 
7 = 1.8 and taking into account the first three roots of 
the corresponding transcendental equations: ji{x) = 
or ni (x)ji (jx) = ni(~fx)ji(x) for i = 1 or 2. In the fig- 
ure we can see the dependence of the radial part of the 
potential on r, and two regions are displayed. For the 
first part (r < a) what we actually have is the electric 
potential $W(r), due to LO phonons of type "1". In 
the second part (a < r < b) the displayed potential is 
$( 2 )(r), due to LO phonons of type "2". In Fig. 1 we 
have used a dotted, a dash-dotted, and a solid line in 
order to distinguish the three first roots of the transcen- 
dental equations in increasing order. As expected, for 
higher values of the corresponding root the electric po- 
tential oscillates more intensively, an effect that should 
led to lower values of the electron-phonon matrix ele- 
ments. On the other hand, it is interesting to note that 
the intensity of the potential, i.e., the absolute value of 



5 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 1.72 1.74 1.76 1.78 1.80 1.82 

radius [units of a] y = b/a 



FIG. 1: Radial part of the electric potential (in units of $o) 
as a function of the radius (in units of a) for I = and 7 = 
1.8. The potentials correspond to the first three LO phonon 
modes and are depicted by a dotted, dash-dotted and solid 
line respectively. 



FIG. 3: Energy levels for the unperturbed CdS/HgS QD/QW 
heterostructure (in eV) as a function of 7 = b/a. For a = 2.5 
nm we show the first two levels (dotted lines). For a — 3.5 
nm (solid lines) the first three energy levels are shown. For 
each level the corresponding value of I is explicitly indicated. 
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FIG. 2: Same as in Fig. 1 but now displaying the potentials 
for the first three LO phonon modes with 1 = 1. 

$«(r), is lower within the core sphere (r < a) in com- 
parison with the layer region a < r < b. When we limit 
ourselves to a discussion of just the radial part of the po- 
tential, we are loosing some information concerning the 
angular dependence that is actually non trivial. However, 
the radial part provides a certain estimation of the poten- 



tial strength that gives an insight into the fundamental 
characteristics of the electron-LO phonon interaction. In 
Fig. 2 we analyze the radial part of the potential taking 
now I — 1 and the same value for 7. Again the first three 
roots of the transcendental equations were considered for 
this value of Z. It is evident from the figure that, in 
contrast with the results shown in Fig. 1, in the present 
case the potential intensity is higher inside the core re- 
gion. From the latter results it could be inferred that the 
electron should feel a stronger phonon interaction in the 
region r < a for 1 = 1, while for I = the electron-phonon 
interaction should be stronger in the region a < r < b. 
However, this kind of analysis could be certainly mis- 
leading. In order to properly understand the strength of 
the electron-phonon interaction we must also take into 
account the electron wave functions. For instance, if for 
I = we get a larger probability amplitude of finding the 
electron inside the core sphere, then this region could 
provide the larger contribution to the electron-phonon 
interaction. For that reason the strength of the inter- 
action is better revealed in such quantities as relaxation 
rates, polaronic corrections, etc., where the matrix el- 
ements of the electron-phonon interaction Hamiltonian 
must be calculated. 

In Fig. 3 we are showing (solid lines) the first three en- 
ergy levels of the QD/QW heterostructure as a function 
of 7 = b/a. These results were obtained from Eq.J33J) 
without the inclusion of the electron-phonon (polaronic) 
corrections. For the core sphere radius we have taken 
a = 3.5 nm, while the band-offset between the conduc- 
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tion bands of CdS and -HgS was estimated as Uq = 1.34 
eV following Ref. [24|. For this value of a the first 
(ground state) energy level is localized inside the QW 
region (E < for our choice of the zero energy level at 
the bottom of the CdS conduction band). The first ex- 
cited energy level (with I — 1) is already out of the QW 
with E > and the second excited level (with I = 0) can 
also be seen in the figure. The values of 7 were chosen in 
a relatively narrow interval 1.72 < 7 < 1.82. According 
to reported experimental result a 1 ! 2 ' 3 ' 4 , 7 is bound to be 
not very large. On the other hand, for lower values of 7 
there are not possible levels for E < 0. For this energy 
level the wave function in the region r < a is expressed 
through hyperbolic functions, while for other energy lev- 
els wavefunctions are given by oscillating trigonometric 
functions. The latter fact provides important changes in 
the calculation of matrix elements and leads to differ- 
ences in polaronic energy shifts (to be discussed below). 
In Fig. 3, using dotted lines, we show the first two en- 
ergy levels for a core sphere radius with a = 2.5 nm. It 
is clear from the figure that, in this latter case, there 
are not energy levels with E < in the analyzed inter- 
val for 7. In Fig. 3 we are not showing the third energy 
level for a = 2.5 nm with I — 0. This level has the energy 
E « 1.30 eV when 7 = 1.72 decreasing down to E « 1.08 
cV for 7 = 1.8. 

As discussed in Section III, the electron-phonon inter- 
action leads to an energy shift of the electronic levels to- 
wards the lower energies, the so-called "polaronic effect" . 
Polaronic effects in the electron, hole, and exciton energy 
levels for QD structures have been studied by different 
authorsiSiiiiiSiSSiSiiSS. In most of the published papers 
the single spherical QD was considered, while the ap- 
plied theoretical formalism is in some cases perturbative 
and in other cases the variational one. For a system like 
the one studied in the current paper, it is plausible to 
apply a perturbative approach as was briefly described 
in Section III. Polaronic corrections may be discussed 
considering both the ground state energy level and the 
excited states levels, and usually involve tedious lengthy 
calculations. In order to make our calculations we ap- 
plied Eq.J22J where the labels p and p' were adapted 
to our system. The electron-LO phonon interaction was 
evaluated according to Ea. (|16fl considering both types 
of LO phonons present in the QD/QW heterostructure 
and including summations over an appropriate amount 
of phonon modes. On the other hand, electron wave- 
functions, properly normalized, were taken according to 
Eq. (|19fl and the discrete energy levels were calculated by 
solving the transcendental equation (|23|l . 

In Fig. 4 we present the absolute value of the pola- 
ronic energy shift corresponding to the case with a = 2.5 
nm. We show the corrections to the ground state en- 
ergy level (the first dotted curve in Fig. 3) for a rela- 
tively wide range of 7 values. These energy shifts, given 
here in meV, should be subtracted from the correspond- 
ing energies reported in Fig. 3. The energy shifts are 
decreasing for growing values of 7, while the calculated 
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y = b/a 

FIG. 4: Absolute value of the energy shifts (in meV) as a 
function of 7 = b/a for the polaronic corrections to the ground 
state (solid line) and excited state (dotted line) energy levels 
considering a QD/QW heterostructure with core radius a — 
2.5 nm; These energy shifts should be subtracted from the 
corresponding energy levels presented in Fig. 3. 

values are relatively low as could be expected for a sys- 
tem having a weak electron-LO phonon coupling. The 
obtained results disagree with those reported in Ref. [2]| 
for a similar situation in which the opposite trend for 
the polaronic corrections was obtained. In Fig. 4 we also 
show the polaronic corrections to the excited energy level 
(dotted curve). This energy level, as mentioned above, 
has I = and was not shown in Fig. 3. The curve, giv- 
ing the dependence with 7, shows a decreasing behavior 
in its first part but, after reaching a minimum approxi- 
mately at 7 = 1.7, shows an increasing behavior for the 
last part of the considered interval. This peculiar kind 
of dependence with the thickness of the HgS layer seems 
to characterize the polaronic corrections for the excited 
levels and was in a certain way commented in Ref. [22| . 
It must also be remarked that, for most of the 7 interval, 
the corrections to the excited level are lower, than those 
to the ground state energy level. 

In Fig. 5 we present the same dependencies as in Fig. 4, 
but taking now a = 3.5. As easily seen in the figure 
the polaronic corrections to the excited state (consider- 
ing again the state with I = 0, i.e., the third solid line 
in Fig. 3) are now considerably higher than the correc- 
tions related to the ground state energy level. Actually, 
the polaronic energy shifts corresponding to the ground 
state energy level arc quite low and difficult to detect 
in the experimental measurements. The latter results, 
in sharp contrast with those described in Fig. 4, can be 
easily understood from the physical point of view. The 
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3- 

2 







excited state energy 
core radius = 3.5 nm 



ground state energy 
core radius = 3.5 nm 



— i — 1 — i — 1 — i — r 

1.77 1.78 1.79 1.80 



1.72 1.73 1.74 1.75 1.76 

y = b/a 



FIG. 5: Same as in Fig. 4 but taking a = 3.5 nm. In both 
cases the results correspond to the CdS/HgS QD/QW het- 
erostructure. 



where Eq.© was applied and /^( r ) are the two func- 
tions given in Eqs. J7J and JSJ. In the above expressions, 
for the sake of brevity, we are avoiding some subscripts. 
The volumes Vi are: for i — 1 the sphere of radius a 
and for i = 2 the spherical shell defined by a < r < b. 
Integrating by parts the integral in Ea. (jAlfl . we first ob- 
tain a surface integral that is always zero because of the 
boundary conditions fulfilled by the mentioned functions 
(notice that these functions are zero over the spherical 
surfaces defining the boundaries of volumes V\ and V2 
respectively). Taking into account that 



2f(») 



the integral can be casted as 



2 f (i) 



-q z f 



f V/ (i) * • V/ (i) d 3 r = q 2 [ /»' 
JVi JVi 



f^d 3 



(A2) 



(A3) 



The latter integral contains an integration over the an- 
gular part of the functions that is easily evaluated in the 
form: J Y* m Yi m dVL — 1. There remain two integrals over 
the radial parts that are different for i = 1, 2 and given 
by 



ground state wave function (for a = 3.5 nm) in the region 
r < a is given by hyperbolic functions (in the framework 
of the used model-system), thus leading to very low prob- 
ability amplitudes if compared with states having E > 0. 
We must also notice that in Fig. 4 we are considering a 
narrower interval for 7 in order to ensure the presence of 
the ground state energy with E < 0. 

As general conclusions to the present work, we can 
emphasize that polaronic corrections are actually inter- 
esting in the study of the electronic energy spectra of 
QD/QW heterostructures. They can be detected in the 
appropriate optical experiments, differences between the 
ground state and excited states polaronic corrections may 
be studied and should provide good estimations of the 
electron-LO phonon coupling in this kind of heterostruc- 
tures. Concerning degenerate energy levels (not studied 
in the current work) the electron-LO phonon interaction 
could introduce detectable energy splittings between the 
degenerate levels, a subject that requires further analysis. 

Appendix A. Normalization of the electric poten- 
tial 

When we substitute the hermitized operator version of 
it 2 , i.e., u 2 — > (tr -u+u-v/)/2, in the classical expression 
for the kinetic energy due to the ion vibrations, we are 
directly led to: 

fig = 1^4, jf V/»* • V/W d\ (at a + 1/2) , 

(Al) 



£jf(qrydr=^jf +1 (qa), (A4) 



/ [ni(qa)ji(qr) - ji{qa)m{qr)f r 2 dr = —Mi(qa) 

(A5) 

where Mi(x) was given in Ea. l|ll|l . For the integral (|A4|) 
the condition ji{qa) = is assumed, while for the inte- 
gral (|A5|I the corresponding condition is ni(qa)jt(qb) = 
ji{qa)ni{qb) . The latter results, together with the analy- 
sis made in Section II, lead to the values of (reported 
in Eqs. (|5J and l|10|l). Using now the relation 

Af\q) = yfAir/HUiL/e^v^ (A6) 

we finally obtain the normalization constants for the elec- 
tric potential reported in Eqs. I|12|l and (|13(l . 
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